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Abstract 

The nonbaryonic dark matter of the Universe is assumed to consist 
of new stable particles. A specific case is possible, when new stable 
particles bear ordinary electric charge and bind in heavy " atoms" by 
ordinary Coulomb interaction. Such possibility is severely restricted 
by the constraints on anomalous isotopes of light elements that form 
positively charged heavy species with ordinary electrons. The trouble 
is avoided, if stable particles X with charge -2 are in excess over 
their antiparticles (with charge +2) and there are no stable particles 
with charges +1 and -1. Then primordial helium, formed in Big Bang 
Nucleosynthesis, captures all X in neutral "atoms" of O-helium 
(OHe). Schrodinger equation for system of nucleus and OHe is con- 
sidered and reduced to an equation of relative motion in a spherically 
symmetrical potential, formed by the Yukawa tail of nuclear scalar 
isoscalar attraction potential, acting on He beyond the nucleus, and 
dipole Coulomb repulsion between the nucleus and OHe at small dis- 
tances between nuclear surfaces of He and nucleus. The values of 
coupling strength and mass of c-meson, mediating scalar isoscalar 
nuclear potential, are rather uncertain. Within these uncertainties 
and in the approximation of rectangular potential wells and wall we 
find a range of these parameters, at which the sodium nuclei have a 
few keV binding energy with OHe. The result also strongly depend 
on the precise value of parameter d that determines the size of nu- 
clei. At nuclear parameters, reproducing DAMA results, OHe-nucleus 
bound states can exist only for intermediate nuclei, thus excluding di- 
rect comparison with these results in detectors, containing very light 
(e.g. 3 He) and heavy nuclei (like Xe). 



1 Introduction 

Ordinary matter around us consists of neutral atoms, in which electrically 
charged nuclei are bound with electrons. Ordinary matter is luminous 
because of electron transitions in atoms. It is stable owing to stability of 
its constituents. Electron is the lightest charged particle. It is stable due 
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to conservation of electromagnetic charge that reflects local gauge U(l) 
invariance. Electromagnetic charge is the source of the corresponding U(l) 
gauge field, electromagnetic field. Nuclei are stable because of stability of 
nucleons. The lightest nucleon - proton - is the lightest baryon and stable 
due to conservation of baryon charge. There is no gauge field related with 
baryon charge. Therefore there are two examples of stable charged particles 
of the ordinary matter: protected by gauge symmetry and protected by 
conserved charge. This excursus in known physics can give us some idea 
on possible constituents of dark atoms, maintaining the dark matter of the 
Universe. 

According to the modern cosmology, the dark matter, corresponding to 
25% of the total cosmological density, is nonbaryonic and consists of new 
stable particles. Such particles (see e.g. [TJ [2j [3] for review and reference) 
should be stable, saturate the measured dark matter density and decouple 
from plasma and radiation at least before the beginning of matter domi- 
nated stage. The easiest way to satisfy these conditions is to involve neutral 
elementary weakly interacting particles. However it is not the only particle 
physics solution for the dark matter problem and more evolved models of 
self-interacting dark matter are possible. In particular, new stable particles 
may possess new U(l) gauge charges and bind by Coulomb-like forces in 
composite dark matter species. Such dark atoms would look nonluminous, 
since they radiate invisible light of U(l) photons. Historically mirror mat- 
ter (see [TJ 2] for review and references) seems to be the first example of 
such a nonluminous atomic dark matter. 

Glashow's tera-helium |5j has offered a new solution for dark atoms of 
dark matter. Tera-U-quarks with electric charge +2/3 formed stable (UUU) 
+2 charged "clusters" that formed with two -1 charged tera-electrons E neu- 
tral [(UUU)EE] tera-helium "atoms" that behaved like Weakly Interacting 
Massive Particles (WIMPs). The main problem for this solution was to 
suppress the abundance of positively charged species bound with ordinary 
electrons, which behave as anomalous isotopes of hydrogen or helium. This 
problem turned to be unresolvable [6] , since the model [5] predicted stable 
tera-electrons E~ with charge -1. As soon as primordial helium is formed 
in the Standard Big Bang Nucleosynthesis (SBBN) it captures all the free 
E~ in positively charged (HeE) + ion, preventing any further suppression 
of positively charged species. Therefore, in order to avoid anomalous iso- 
topes overproduction, stable particles with charge -1 (and corresponding 
antiparticles) should be absent, so that stable negatively charged particles 
should have charge -2 only. 

Elementary particle frames for heavy stable -2 charged species are pro- 
vided by: (a) stable " antibaryons" UUU formed by anti-£7 quark of fourth 
generation [Zl El El ED] (b) AC-leptons [lQl EU E2], predicted in the exten- 
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sion [TT] of standard model, based on the approach of almost-commutative 
geometry |13j . (c) Technileptons and anti-technibaryons [14] in the frame- 
work of walking technicolor models (WTC) [15]. (d) Finally, stable charged 
clusters of (anti)quarks U5 of 5th family can follow from the ap- 

proach, unifying spins and charges [16]. Since all these models also predict 
corresponding +2 charge antiparticles, cosmological scenario should pro- 
vide mechanism of their suppression, what can naturally take place in the 
asymmetric case, corresponding to excess of -2 charge species, X~~ . Then 
their positively charged antiparticles can effectively annihilate in the early 
Universe. 

If new stable species belong to non-trivial representations of electroweak 
SU(2) group, sphaleron transitions at high temperatures can provide the re- 
lationship between baryon asymmetry and excess of -2 charge stable species, 
as it was demonstrated in the case of WTC [Ti ] H71 H51 [TO] ■ 

After it is formed in the Standard Big Bang Nucleosynthesis (SBBN), 
4 iie screens the X~~ charged particles in composite ( 4 He ++ X~~) O- 
helium "atoms" [8j. For different models of X these "atoms" are also 
called ANO-helium P [TO], Ole-helium [TOJ [TO] or techni-O-helium [TO] . 
We'll call them all O-helium {OHe) in our further discussion, which fol- 
lows the guidelines of [20] . 

In all these forms of O-helium, X~~ behaves either as lepton or as 
specific "heavy quark cluster" with strongly suppressed hadronic interac- 
tion. Therefore O-helium interaction with matter is determined by nuclear 
interaction of He. These neutral primordial nuclear interacting objects con- 
tribute to the modern dark matter density and play the role of a nontrivial 
form of strongly interacting dark matter |21[ [22] . 

Here after a brief review of the qualitative picture of OHe cosmological 
evolution [HI [TOJ [TOJ [TOJ [23] we concentrate on some open questions in the 
properties of these dark atoms and their interaction with matter. This 
analysis is used in our second contribution to explain the puzzles of dark 
matter searches [24] 

2 Some features of O-helium Universe 

Following [HJ [9[ [TOJ [TOJ [TOJ [TOJ [20] consider charge asymmetric case, when 
excess of X provides effective suppression of positively charged species. 
In the period 100 s < t < 300 s at 100 keV > T > T Q = I Q /27 ra 60keV, 
He has already been formed in the SBBN and virtually all free X are 
trapped by He in O-helium "atoms" (^He ++ X ). Here the O-helium 
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ionization potential i^] 

I = ZlZ 2 He a 2 m He /2 ps 1.6 MeV, (1) 

where a is the fine structure constant, Zu e = 2 and Z x = 2 stands for the 
absolute value of electric charge of X . The size of these "atoms" is [8l[T2] 

Ro ~ l/{Z x Z He am He ) « 2 • 1(T 13 cm (2) 

Here and further, if not specified otherwise, we use the system of units 
H = c = k = 1. 

Due to nuclear interactions of its helium constituent with nuclei in the 
cosmic plasma, the O-helium gas is in thermal equilibrium with plasma and 
radiation on the Radiation Dominance (RD) stage, while the energy and 
momentum transfer from plasma is effective. The radiation pressure acting 
on the plasma is then transferred to density fluctuations of the O-helium 
gas and transforms them in acoustic waves at scales up to the size of the 
horizon. 

At temperature T < T od ~ 2005g //3 eV the energy and momentum 
transfer from baryons to O-helium is not effective [8j [T3] because 

n B {crv) (m p /m )t < 1, 

where m is the mass of the OHe atom and 53 = m J (1 TeV). Here 

a ps a ~ itR 2 ps 1(T 25 cm 2 , (3) 



and v = y / 2T '/m p is the baryon thermal velocity. Then O-helium gas 
decouples from plasma. It starts to dominate in the Universe after t ~ 10 12 s 
at T < Trm ~ 1 eV and O-helium "atoms" play the main dynamical role 
in the development of gravitational instability, triggering the large scale 
structure formation. The composite nature of O-helium determines the 
specifics of the corresponding dark matter scenario. 

At T > Trm the total mass of the OHe gas with density pd = (Trm /T)p to t 
is equal to 

4vr 3 4tt T rm mpi 2 
M = —p d t = - — —m Pl { — ) 

within the cosmological horizon 1^ = t. In the period of decoupling T = T 0( i, 
this mass depends strongly on the O-helium mass S3 and is given by |14] 

M od = ^m Pl (^f - 2 • 10 44 5 3 - 2 g = lO^Mo, (4) 

J-od J-od 



1 Thc account for charge distribution in He nucleus leads to smaller value I w 1.3 MeV 
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where Mq is the solar mass. O-helium is formed only at T Q and its total 
mass within the cosmological horizon in the period of its creation is M Q = 
M od {T od /T f = 10 37 g. 

On the RD stage before decoupling, the Jeans length Aj of the OHe 
gas was restricted from below by the propagation of sound waves in plasma 
with a relativistic equation of state p = e/3, being of the order of the 
cosmological horizon and equal to Xj = Ih/VS = tj\f?>. After decoupling 
at T = T od , it falls down to Aj ~ v Q t, where v Q = y / 2T od /m . Though after 
decoupling the Jeans mass in the OHe gas correspondingly falls down 

Mj ~ v 3 Q M od ~ 3 • W~ 14 M od , 

one should expect a strong suppression of fluctuations on scales M < M , as 
well as adiabatic damping of sound waves in the RD plasma for scales M Q < 
M < M od . It can provide some suppression of small scale structure in the 
considered model for all reasonable masses of O-helium. The significance 
of this suppression and its effect on the structure formation needs a special 
study in detailed numerical simulations. In any case, it can not be as strong 
as the free streaming suppression in ordinary Warm Dark Matter (WDM) 
scenarios, but one can expect that qualitatively we deal with Warmer Than 
Cold Dark Matter model. 

Being decoupled from baryonic matter, the OHe gas does not follow 
the formation of baryonic astrophysical objects (stars, planets, molecular 
clouds...) and forms dark matter halos of galaxies. It can be easily seen 
that O-helium gas is collisionless for its number density, saturating galactic 
dark matter. Taking the average density of baryonic matter one can also 
find that the Galaxy as a whole is transparent for O-helium in spite of its 
nuclear interaction. Only individual baryonic objects like stars and planets 
are opaque for it. 

3 Signatures of O-helium dark matter in the Galaxy 

The composite nature of O-helium dark matter results in a number of 
observable effects, which we briefly discuss following [18]. 

3.1 Anomalous component of cosmic rays 

O-helium atoms can be destroyed in astrophysical processes, giving rise to 
acceleration of free X in the Galaxy. 

O-helium can be ionized due to nuclear interaction with cosmic rays 
[8j [20]. Estimations [8j [25] show that for the number density of cosmic 
rays ucr = 10 -9 cm -3 during the age of Galaxy a fraction of about 10~ 6 
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of total amount of OHe is disrupted irreversibly, since the inverse effect of 
recombination of free X is negligible. Near the Solar system it leads to 
concentration of free X nx = 3 • 10 _10 5^~ 1 cm -3 . After OHe destruction 

-1 Icy 

free X~~ have momentum of order px — \/2 • Mx • I — 2GeVS 3 and 
velocity v/c = 2 ■ lO -3 ^ and due to effect of Solar modulation these 
particles initially can hardly reach Earth [17\ 125] . Their acceleration by 
Fermi mechanism or by the collective acceleration forms power spectrum 
of X component at the level of X/p ~ nx/n g = 3 • lO -10 ^ 1 , where 
n g ~ 1 cm~ 3 is the density of baryonic matter gas. 

At the stage of red supergiant stars have the size ~ 10 15 cm and during 
the period of this stage~ 3- 10 15 s, up to ~ lO -9 ^ 1 of O-helium atoms per 
nucleon can be captured [T71 [25] . In the Supernova explosion these OHe 
atoms are disrupted in collisions with particles in the front of shock wave 
and acceleration of free X by regular mechanism gives the corresponding 
fraction in cosmic rays. However, this picture needs detailed analysis, based 
on the development of OHe nuclear physics and numerical studies of OHe 
evolution in the stellar matter. 

If these mechanisms of X acceleration are effective, the anomalous 
low Z/A component of —2 charged X can be present in cosmic rays at 
the level X/p ~ n x /n g ~ lCT 9 ^ 1 , and be within the reach for PAMELA 
and AMS02 cosmic ray experiments. 

In the framework of Walking Technicolor model the excess of both stable 
X~~ and Y ++ is possible [17], the latter being two-three orders of mag- 
nitude smaller, than the former. It leads to the two-component composite 
dark matter scenario with the dominant OHe accompanied by a subdomi- 
nant WIMP-like component of (X~~Y ++ ) bound systems. Technibaryons 
and technileptons can be metastable and decays of X and Y ++ can pro- 
vide explanation for anomalies, observed in high energy cosmic positron 
spectrum by PAMELA and in high energy electron spectrum by FERMI 
and ATIC. 

3.2 Positron annihilation and gamma lines in galactic bulge 

Inelastic interaction of O-helium with the matter in the interstellar space 
and its de-excitation can give rise to radiation in the range from few 
keV to few MeV. In the galactic bulge with radius ~ 1 kpc the num- 
ber density of O-helium can reach the value n Q ~ 3 • 1CP 3 /S3 cm -3 and 
the collision rate of O-helium in this central region was estimated in [20J: 
dN/dt = nlavhAnrl/Z ps 3-10 42 5^ 2 s -1 . At the velocity of v h ~ 3-10 7 cm/s 
energy transfer in such collisions is AE ~ 1 MeV 53. These collisions can 
lead to excitation of O-helium. If 2S level is excited, pair production dom- 
inates over two-photon channel in the de-excitation by EO transition and 
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positron production with the rate 3 • 10 42 5 3 ^ 2 s~ 1 is not accompanied by 
strong gamma signal. According to [26J this rate of positron production for 
S3 ~ 1 is sufficient to explain the excess in positron annihilation line from 
bulge, measured by INTEGRAL (see [27] for review and references). HOHe 
levels with nonzero orbital momentum are excited, gamma lines should be 
observed from transitions (n > m) E nm = 1.598 MeV(l/m 2 — 1/n 2 ) (or 
from the similar transitions corresponding to the case I a = 1.287 MeV) at 
the level 3 • W~ 4 S^ 2 ( cm 2 sMeVster)" 1 . 

It should be noted that the nuclear cross section of the O-helium interac- 
tion with matter escapes the severe constraints [22] on strongly interacting 
dark matter particles (SIMPs) |21 j. [22] imposed by the XQC experiment 
[28] , Therefore, a special strategy of direct O-helium search is needed, as 
it was proposed in [29] . 



4 O-helium interaction with nuclei 

The evident consequence of the O-helium dark matter is its inevitable pres- 
ence in the terrestrial matter, which appears opaque to O-helium and stores 
all its in-falling flux. After they fall down terrestrial surface, the in-falling 
OHe particles are effectively slowed down due to elastic collisions with mat- 
ter. In underground detectors, OHe "atoms" are slowed down to thermal 
energies and give rise to energy transfer ~ 2.5 • 10 -4 eVA/S^, far below the 
threshold for direct dark matter detection. It makes this form of dark mat- 
ter insensitive to the severe CDMS constraints [30j. However, OHe induced 
processes in the matter of underground detectors can result in observable 
effects. These effects, considered in a separate contribution [M], strongly 
depend on the details of the OHe interaction with nuclei, which we consider 
here. 



4.1 Structure of X atoms with nuclei 

The properties of OHe interaction with matter are determined first of all 
by the structure of OHe atom that follows from the general analysis of the 
bound states of X with nuclei. 

Consider a simple model |31j . in which the nucleus is regarded as a 
sphere with uniform charge density and in which the mass of the X 
is assumed to be much larger than that of the nucleus. Spin dependence 
is also not taken into account so that both the particle and nucleus are 
considered as scalars. Then the Hamiltonian is given by 

p 2 ZZ x a ZZ x a . r 2 
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for short distances r < R and 



H = -g--?%°, (6) 
2Am p R ' w 

for long distances r > R, where a is the fine structure constant, R = 
doA 1 / 3 ~ 1.2A 1 / 3 / (20QMeV) is the nuclear radius, Z is the electric charge 
of nucleus and Z x = 2 is the electric charge of negatively charged particle 
X . Since Am p <C Mx the reduced mass is 1/m = l/(Am p ) + 1/Mx ~ 
l/(Am p ). 

For small nuclei the Coulomb binding energy is like in hydrogen atom 
and is given by 

E b = ^Z 2 Z 2 x a 2 Am p . (7) 

For large nuclei X is inside nuclear radius and the harmonic oscillator 
approximation is valid for the estimation of the binding energy 

_ 2>,ZZ x a 1 , ZZ x a . i /9 . 

For the intermediate regions between these two cases with the use of 
trial function of the form if) ~ eT^ r l R variational treatment of the problem 
[3T] gives 

Eb= A^m F{ZZxaAmpRl (9) 
where the function F{a) has limits 

F(a -> 0) -> -a 2 - -a 4 (10) 
2 5 

and 

F(a -»• oo) -)■ - (3a) 1/2 , (11) 

where a = Z Z x aAm p R. For < a < 1 the Coulomb model gives a good 
approximation, while at 2 < a < oo the harmonic oscillator approximation 
is appropriate. 

In the case of OHe a = ZZ x aAm p R < 1, what proves its Bohr-atom-like 
structure, assumed in our earlier papers [U [101 HH [181 H3 [20]. However, 
the size of He, rotating around X in this Bohr atom, turns out to be of 
the order and even a bit larger than the radius r Q of its Bohr orbit, and 
the corresponding correction to the binding energy due to non-point-like 
charge distribution in He is significant. 

Bohr atom like structure of OHe seems to provide a possibility to use 
the results of atomic physics for description of OHe interaction with matter. 
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However, the situation is much more complicated. OHe atom is similar to 
the hydrogen, in which electron is hundreds times heavier, than proton, so 
that it is proton shell that surrounds "electron nucleus". Nuclei that in- 
teract with such "hydrogen" would interact first with strongly interacting 
"protonic" shell and such interaction can hardly be treated in the frame- 
work of perturbation theory. Moreover in the description of OHe interaction 
the account for the finite size of He, which is even larger than the radius 
of Bohr orbit, is important. One should consider, therefore, the analysis, 
presented below, as only a first step approaching true nuclear physics of 
OHe. 



4.2 Potential of O-helium interaction with nuclei 

Our explanation PIQUES] of the results of DAMA/Nal [33] and DAM A/LIBRA 
|34j experiments is based on the idea that OHe, slowed down in the mat- 
ter of detector, can form a few keV bound state with nucleus, in which 
OHe is situated beyond the nucleus. Therefore the positive result of these 
experiments is explained by reaction 

A + ( 4 Fe ++ A— ) [A( 4 He ++ X—)] + 7 (12) 

with nuclei in DAMA detector. 

In our earlier studies [HI [191 [32] the conditions were found, under which 
both sodium and iodine nuclei have a few keV bound states with OHe, 
explaining the results of DAMA experiments by OHe radiative capture to 
these levels. Here we extend the set of our solutions by the case, when the 
results of DAMA experiment can be explained by radiative OHe capture 
by sodium only and there are no such bound states with iodine and Tl. 

Schrodinger equation for OHe-nucleus system is reduced (taking apart 
the equation for the center of mass) to the equation of relative motion for 
the reduced mass 

Am p m 

m = — , (13) 

Atrip + m 

where m p is the mass of proton and m ~ Mx + 4m p is the mass of 
OHe. Since m ~ Mx 3> Am p , center of mass of Ohe-nucleus system 
approximately coincides with the position of X . 

In the case of orbital momentum /=0 the wave functions depend only 
on r. 

The approach of [181 021 E2] assumes the following picture: at the dis- 
tances larger, than its size, OHe is neutral, being only the source of a 
Coulomb field of X screened by He shell 

U c = Z * Za - F ^ r \ (14) 
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where Zx = — 2 is the charge of X , Z is charge of nucleus, Fx(r) = 
(1 + r /r )exp{— 2r/r Q ) is the screening factor of Coulomb potential (see 
e.g. [35]) of X and r Q is the size of OHe. Owing to the negative sign of 
Zx = —2, this potential provides attraction of nucleus to OHe. 

Then helium shell of OHe starts to feel Yukawa exponential tail of at- 
traction of nucleus to He due to scalar-isoscalar nuclear potential. It should 
be noted that scalar-isoscalar nature of He nucleus excludes its nuclear in- 
teraction due to 7r or p meson exchange, so that the main role in its nuclear 
interaction outside the nucleus plays a meson exchange, on which nuclear 
physics data are not very definite. The nuclear potential depends on the 
relative distance between He and nucleus and we take it in the form 

jj = A He Ag 2 exp(-p\r- p\) 
\f-p\ 

Here r is radius vector to nucleus, p is the radius vector to He in OHe, 
Af{ e = 4 is atomic weight of helium, A is atomic weight of nucleus, p and 
g 2 are the mass and coupling of a meson - mediator of nuclear attraction. 

Strictly speaking, starting from this point we should deal with a three- 
body problem for the system of He, nucleus and X~~ and the correct 
quantum mechanical description should be based on the cylindrical and 
not spherical symmetry. In the present work we use the approximation of 
spherical symmetry and take into account nuclear attraction beyond the 
nucleus in a two different ways: 1) nuclear attraction doesn't influence the 
structure of OHe, so that the Yukawa potential (fT5j) is averaged over \p\ 
for spherically symmetric wave function of He shell in OHe; 2) nuclear 
attraction changes the structure of OHe so that He takes the position \p\ = 
r Q , which is most close to the nucleus. Due to strong attraction of He by 
the nucleus the second case (which is labeled "b" in successive numerical 
calculations) seems more probable. In the lack of the exact solution of the 
problem we present both the results, corresponding to the first case (which 
are labeled "m" in successive numerical calculations), and to the second 
case (which is labeled "b") in order to demonstrate high sensitivity of the 
numerical results to choice of parameters. 

In the both cases nuclear attraction results in the polarization of OHe 
and the mutual attraction of nucleus and OHe is changed by Coulomb 
repulsion of He shell. Taking into account Coulomb attraction of nucleus 
by X one obtains dipole Coulomb barrier of the form 

Z He Zar 

U d = ~2 • (16) 

When helium is completely merged with the nucleus the interaction 
is reduced to the oscillatory potential JI7J) of X with homogeneously 
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Figure 1: The approximation of rectangular well for potential of OHe- 
nucleus system. 



charged merged nucleus with the charge Z + 2, given by 

3, (Z + 2)Z x a 1 ]Z + 2)^q 1/2 
^ m = 2 ( 1 R^A + A^R* } " (17) 

To simplify the solution of Schrodinger equation we approximate the 
potentials (fT4"|) - (fTT|) by a rectangular potential that consists of a potential 
well with the depth U\ at r < c = R, where R is the radius of nucleus, of a 
rectangular dipole Coulomb potential barrier t/ 2 at R<r<a = R + r + 
The, where r\ ie is radius of helium nucleus, and of the outer potential well 
U3, formed by the Yukawa nuclear interaction (|15p and residual Coulomb 
interaction (1141) , The values of U\ and U2 were obtained by the averaging 
of the (|17p and (JTSJ) in the corresponding regions, while U3 was equal to 
the value of the nuclear potential (|15p at r = a and the width of this outer 
rectangular well (position of the point b) was obtained by the integral of the 
sum of potentials (|15p and (|14p from a to 00. It leads to the approximate 
potential, presented on Fig. [H 

Solutions of Schrodinger equation for each of the four regions, indicated 
on Fig. [H are given in textbooks (see e.g.[35j) and their sewing determines 
the condition, under which a low-energy OHe-nucleus bound state appears 
in the region III. 
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4.3 Low energy bound state of O-helium with nuclei 

The energy of this bound state and its existence strongly depend on the 
parameters p and g 2 of nuclear potential (|15p . On the Fig. [2] the regions 
of these parameters, giving 4 keV energy level in OHe bound state with 
sodium are presented. Radiative capture to this level can explain results 
of DAMA/Nal and DAMA/LIBRA experiments with the account for their 
energy resolution [36]. The lower shaded region on Fig. [2] corresponds to 
the case of nuclear Yukawa potential U^ m , averaged over the orbit of He 
in OHe, while the upper region corresponds to the case of nuclear Yukawa 
potential U^b with the position of He most close to the nucleus at p = r D .The 
result is also sensitive to the precise value of d Q , which determines the size 
of nuclei R = d A 1//3 . The two narrow strips in each region correspond 
to the experimentally most probable value d Q = 1.2/(200 MeV). In these 
calculations the mass of OHe was taken equal to m = ITeV, however the 
results weakly depend on the value of m > ITeV. 

4.4 Energy levels in other nuclei 

The important qualitative feature of the presented solution is the restricted 
range of intermediate nuclei, in which the OHe-nucleus state beyond nuclei 
is possible. For the chosen range of nuclear parameters, reproducing the 
results of DAMA/Nal and DAMA/LIBRA, we can calculate the binding 
energy of OHe-nucleus states in nuclei, corresponding to chemical compo- 
sition of set-ups in other experiments. It turns out that there are no such 
states for light and heavy nuclei. In the case of nuclear Yukawa potential 
U^bi corresponding to the position of He most close to the nucleus at p = r , 
the range of nuclei with bound states with OHe corresponds to the part 
of periodic table between B and Ti. This result is stable independent on 
the used scheme of numerical calculations. The upper limits on the nuclear 
parameters p and g 2 , at which there exists OHe-nucleus bound state are 
presented for this case on Figj3j In the case of potential Uz m , averaged 
over the orbit of He in OHe, there are no OHe bound states with nuclei, 
lighter than Be and heavier than Ti. However, the results are very sensitive 
to the numerical factors of calculations and the existence of OHe-Ge and 
OHe-Ga bound states at a narrow window of parameters p and g 2 turns to 
be strongly dependent on these factors so that change in numbers smaller 
than 1% can give qualitatively different result for Ge and Ga. The results 
for the case (m) are shown on FigJH Both for the cases (b) and (m) there 
is a stable conclusion that there are no OHe-nucleus bound states with Xe, 
I and TI. 

For the experimentally preferred value d Q = 1.2/(200 MeV) the results 
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Figure 2: The region of parameters // and g 2 , for which Na-OHe system has 
a level in the interval 4 keV. Two lines determine at d Q = 1.2/(200 MeV) 
the region of parameters, at which the bound system of this element 
with OHe has a 4 keV level. In the region between the two strips the 
energy of level is below 4 keV. There are also indicated the range of 
g 2 /^ 2 (dashed lines) as well as their preferred values (thin lines) de- 
termined in [37] from parametrization of the relativistic (a — u) model 
for nuclear matter. The uncertainty in the determination of parameter 
1.15/(200 MeV) < d Q < 1.3/(200 MeV) results in the uncertainty of (j, and 
g 2 shown by the shaded regions surrounding the lines. The case of nuclear 
Yukawa potential U% m , averaged over the orbit of He in OHe, corresponds 
to the lower lines and shaded region, while the upper lines and shaded re- 
gion around them illustrate the case of nuclear Yukawa potential U^b with 
the position of He most close to the nucleus at p = r a . 
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Figure 3: Existence of low energy bound states in OHe-nucleus system in 
the case b for nuclear Yukawa potential U^b with the position of He most 
close to the nucleus at p = r a . The lines, corresponding to different nuclei, 
show the upper limit for nuclear physics parameters fj, and g 2 , at which 
these bound states are possible. The choice of parameters corresponding 
to 4 keV OHe-Na bound state, excludes region below Na line. 
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Figure 4: Existence of low energy bound states in OHe-nucleus system in 
the case m for nuclear Yukawa potential U^ m , averaged over the orbit of He 
in OHe. The lines, corresponding to different nuclei, show the upper limit 
for nuclear physics parameters fj, and g 2 , at which these bound states are 
possible. The choice of parameters corresponding to 4 keV OHe-Na bound 
state, excludes region below Na line. 
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Figure 5: Energy levels in OHe bound system with carbon, oxygen, fluorine, 
argon, silicon, aluminium and chlorine for the case of the nuclear Yukawa 
potential U^. The predictions are given for the range of g 2 /^ 2 determined 
in [37] from parametrization of the relativistic (cr — oS) model for nuclear 
matter. The preferred values of g 2 / 'fi 2 are indicated by the corresponding 
marks (squares or circles) 



of calculation of the binding energy of OHe-nucleus systems for carbon, 
oxygen, fluorine, argon, silicon, aluminium and chlorine are presented on 
Fig. [5] for the case of the nuclear Yukawa potential U^b and on Fig. [6] for 
the case of the potential U^ m . The difference in these results demonstrates 
their high sensitivity to the choice of parameters. 



5 Conclusions 

To conclude, the results of dark matter search in experiments DAMA/Nal 
and DAMA/LIBRA can be explained in the framework of composite dark 
matter scenario without contradiction with negative results of other groups. 
This scenario can be realized in different frameworks, in particular in Min- 
imal Walking Technicolor model or in the approach unifying spin and 
charges and contains distinct features, by which the present explanation 
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Figure 6: Energy levels in OHe bound system with carbon, oxygen, fluorine, 
argon, silicon, aluminium and chlorine for the case of the nuclear Yukawa 
potential U^ m . The predictions are given for the range oi g 2 / fi 2 determined 
in [37] from parametrization of the relativistic (a — oS) model for nuclear 
matter. The preferred values of g 2 / fi 2 are indicated by the corresponding 
marks (squares or circles) 
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can be distinguished from other recent approaches to this problem [38] (see 
also review and more references in |39|). 

Our explanation is based on the mechanism of low energy binding of 
OHe with nuclei. We have found that within the uncertainty of nuclear 
physics parameters there exists their range at which OHe binding energy 
with sodium is equal to 4 keV and there is no such binding with iodine and 
thallium. 

With the account for high sensitivity of our results to the values of 
uncertain nuclear parameters and for the approximations, made in our cal- 
culations, the presented results can be considered only as an illustration of 
the possibility to explain effects in underground detectors by OHe binding 
with intermediate nuclei. However, even at the present level of our studies 
we can make a conclusion that effects of such binding should strongly differ 
in detectors with the content, different from Nal, and can be absent in de- 
tectors with very light (e.g. 3 He) and heavy nuclei (like xenon and probably 
germanium). Therefore test of results of DAMA/Nal and DAMA/LIBRA 
experiments by other experimental groups can become a very nontrivial 
task. 
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